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SUMMARY 


Solutions  are  presented  for  the  lift  and  pitching  moment  on  a  two- 
dimensional  flat  plate  undergoing  a  step  variation  in  pitch  or  vertical 
displacement.  Resulting  indicial  lift  and  pitching-moment  curves  are 
given  for  free-stream  Mach  numbers  equal  to  0,  0.8,  1.0,  1.2,  and  2.0. 
Considerable  use  is  made  of  the  analogy  between  the  boundary  values  for 
a  two-dimensional  wing  in  unsteady  motion  and  those  for  a  three- 
dimensional,  lifting  surface  in  steady  motion.  The  incompressible, 
unsteady  case,  for  which  Wagner*s  classical  treatment  already  exists,  is 
shown  to  be  analogous  to  a  problem  in  slender  wing  theory. 


INTRODUCTION 


The  response  in  lift  for  a  two-dimensional  airfoil  that  starts  sud¬ 
denly  from  rest  and  moves  forward  at  a  constant  velocity  in  an  incom¬ 
pressible  fluid  was  studied  originally  by  Wagner  (reference  1).  It  is 
well  known  that  such  a  response,  referred  to  as  the  indicial  lift  due  to 
angle  of  attack,  can,  by  suitable  superposition,  be  used  to  find  the  lift 
on  a  wing  undergoing  an  arbitrary  variation  of  angle  of  attack  with  time. 
In  particular,  Garrick  (reference  2)  showed  that  Wagner fs  results  were 
consistent  with  those  obtained  by  Theodorsen  (reference  3)  for  a  harmoni¬ 
cally  oscillating  wing. 

Subsequently,  aerodynami cists  have  become  interested  in  the  effect 
of  compressibility  on  the  unsteady  motion  of  a  wing.  Most,  of  this  inter¬ 
est,  however,  has  been  focused  on  the  problem  of  wing  flutter  which 
involves  harmonic  motion.  Hence,  the  analysis  of  the  compressibility 
effects  has  usually  commenced  with  the  immediate  assumption  that  the 
motion  is  harmonic.  Such  studies  which,  in  general,  can  be  considered  as 
an  extension  of  the  work  of  Theodorsen,  have  been  particularly  successful 
in  the  case  of  supersonic  flow  (e.g. ,  reference  4).  In  subsonic  studies, 
however,  these  methods  have  proven  to  be  somewhat  cumbersome. 
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.The  purpose  of  this  report  is  to  calculate  the  effects  of  compress¬ 
ibility  on  the  indicial  responses.  At  supersonic  speeds  (e.g.,  refer¬ 
ence  5)  some  of  these  indicial  curves  have  been  obtained.  The  material 
presented  here,  however,  vill  provide  a  complete  set  of  indicial 
responses  (c^,  c^*,  cz^f,  c^*)  for  supersonic  speeds  and,  what  is 

more  important,  vill  also  provide  these  same  responses  for  subsonic 
speeds;  the  particular  case  of  a  free— stream  Mach  number  equal  to  0.8 
having  been  chosen  for  detailed  consideration.  Results  are  also  pre¬ 
sented  for  a  wing  traveling  at  the  speed  of  sound. 

The  analysis  required  to  formulate  these  indicial  curves  is  based 
on  the  two-dimensional  vave  equation  having  boundary  conditions  speci¬ 
fied  over  the  plane  in  which  the  wing  moves.  The  particular  boundary- 
value  problems  involved  are,  in  fact,  completely  analogous  to  certain 
three-dimensional,  steady-state,  supersonic,  lifting-surface  problems. 

The  latter  field  has  received  much  attention  during  the  last  few  years 
and  analytical  techniques  have  been  developed  by  means  of  which  many 
supersonic  lifting-surface  problems  have  been  solved.  Because  of  the 
anaiogy,  these  techniques  can  be  applied  to  the  field  of  two-dimensional, 
unsteady  flow.  The  analogy  even  extends  to  the  point  where  a  two- 
dimensional  wing  moving  unsteadily  in  an  incompressible  field  has  for 
its  analog  a  three-dimensional,  lifting  surface  moving  steadily  in  a 
compressible  medium  at  the  speed  of  sound.  By  means  of  this  correspond¬ 
ence,  Wagner ? s  classical  result  can  be  rederived  as  a  problem  in  slender- 
ving  theory. 

The  part  of  the  analysis  pertaining  to  the  response  of  an  unsteady 
wing  traveling  at  a  subsonic  speed  is  lengthy  and  somewhat  tedious 
regardless  of  the  method  of  approach.  With  the  use  of  indicial  functions , 
however,  the  calculations  are  reasonably  straightforward,  especially  for 
Mach  numbers  around  0.8  to  1.0.  Further,  the  use  of  indicial  functions 
sheds  considerable  light  on  the  manner  in  which  Mach  number  variations 
affect  the  section  aerodynamic  characteristics. 


LIST  OF  IMPORTANT  SYMBOLS 


m 


'CL 


velocity  of  sound  in  the  free  stream 
chord  of  wing 

section  lift  coefficient 


(*) 
(=0 


section  moment  coefficient 


indicial  section  lift  coefficient  due  to  angle-of-attack 
change  (without  pitching) 


A 


r 


A 
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%' 

indicial  section  lift  coefficient  due  to  pitching  on  a  wing 
rotating  about  its  leading  edge 

r*  » 

C]% 

indicial  section  pitching-moment  coefficient  due  to  angle— of— 
attack  change  (without  pitching)  measured  about  the  leading 
edge  and  considered  positive  when  the  trailing  edge  is 
forced  down 

V 

indicial  section  pitching-moment  coefficient 
measured  about  the  leading  edge  for  a  wing 
its  leading  edge;  considered  positive  when 
edge  is  forced  down 

due  to  pitching 
rotating  about 
the  trailing 

(V) 

3/4 

indicial  section  pitching-moment  coefficient 
measured  about  the  leading  edge  for  a  wing 
its  3 /4-chord  point 

due  to  pitching 
rotating  about 

2 

section  lift 

m 

section  moment 

Mo 

flight  Mach  number 

P 

static  pressure 

Ap 

q. 

P2  “  Pu  : 

6c 

Vo 

dynamic  pressure  (  ip0V02J  * 

r,s 

rd 

oblique  coordinates  f  r  =  s  =  \ 

\  J 

s(l+M0)  -  cj2 

1-Mo 

t» 

time 

t 

a0t« 

t 

o 

a0t*  _  t 
c  c 

u,v,w 

perturbation  velocity  components  in  the  x,y 

yz  directions. 

respectively 


Au 


“  u2 
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VQ  flight  velocity 

x,y,z  Cartesian  coordinates 

x  2. 

X°  c 

a  angle  of  attack  in  radians 

*  J\  ^  I 

r  circulation  strength 

5  function,  normalized  with  respect  to  t1,  thus 

S ( t 1 )dt 1  =  1 

6(t») 

of  pitch  in  radians 
pQ  density  in  undisturbed  air 

cp  perturbation  velocity  potential 

Subscripts 


u 

upper 

side 

of 

z=0 

plane 

? 

lower 

side 

of 

z=0 

plane 

i  variable  of  integration 

5  part  of  a  response  multiplied  by  a  6  function 


6(t») 


Dirac 


6(o) 

e 


lim 
t  f->0 

angle 


METHOD 


A  basic  linearized  form  of  the  partial  differential  equation  which 
governs  the  flow  field  surrounding  a  thin  wing  moving  through  the  air 
can  be  written  in  terms  of  the  perturbation  velocity  potential,  cp ,  as 


^XX  +  ^inr  +  CP- 


yy 


zz 


aQ2 


(i) 
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where 

x,y,z 

spatial  coordinates 

t* 

true  time 

ao 

speed  of  sound 

Equation  (l)  is  applicable  when  the  fluid  at  infinity  is  at  rest  with 
respect  to  the  x,y,z  coordinate  system  and  the  wing  or  body  traces 
certain  space  curves  for  which  time  is  the  parameter;  moreover  the  per¬ 
turbation  velocity  components  u=cpx,  v=cfy,  w=cpz  must  be  small  relative 
to  the  flight  velocity  VQ  of  the  wing. 

In  the  following  analysis  a  flat— plate  wing  moving  in  the  z=0 
plane  away  from  the  origin  along  the  negative  x  axis  will  be  studied. 
Since  the  variations  in  the  flow  along  the  span  will  be  neglected  and 
the  slope  of  the  stream  lines  at  the  wing  surface  will  be  taken  as  the 
ratio  of  wu,  the  vertical  induced  velocity  in  the  z=0  plane,1  to  VQ, 
the  constant  forward  velocity  of  the  wing,  the  analysis  is  consistent 
with  the  assumptions  of  two-dimensional  thin-airfoil  theory.  Hence,  for 
the  incompressible  case,  the  partial  differential  equation  reduces  to 

<Pxx  +  <Pzz  =  0  (2) 


since  the  speed  of  sound  is  very  large  with  respect  to  the  velocity  of 
the  wing,  and  for  the  compressible  case  it  reduces  to 

V  +  <tzz  -’’tt  (3) 

where  t  equals  a0t'. 

The  boundary  conditions  reduce  to  the  specification  of  wu  over 
that  portion  of  the  xy  plane  occupied  by  the  wing  as  time  progresses. 


The  solution  in  terms  of  cp  can  be  converted  into  the  loading 
coefficient  by  means  of  the  equation 

Ap  _  _  k  Sep 

lo  V02  St'  VqMq  St 


00 


The  results  for  section  lift  and  pitching  moment  can  then  be  evaluated 
by  appropriate  integrations. 

■•The  subscript  u  indicates  the  value  of  w  as  z  approaches  zero 
from  above.  Since  the  discussion  is  limited  to  flat— plate  surfaces 
without  thickness ,  this  is  the  same  as  the  value  of  w  obtained 
as  z  approaches  zero  from  below,  and  thus  wu  simply  refers  to  the 
value  of  w  in  the  z=0  plane. 
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Boundary  Conditions  for 


(a) 

namely: 


the  Indicial  Functions 

i 

The  fundamental  boundary— 
value  problem  to  be  considered  in 
this  report  is  the  one  generating 
the  so-called  indicial  curves  for 
loading,  lift,  and  pitching  moment. 
By  definition,  an  indicial  function 
is  the  response  to  a  disturbance 
which  is  applied  abruptly  at  time 
equals  zero  and  is  held  constant 
thereafter;  that  is,  a  disturbance 
given  by  a  step  function.  For 
example,  if  the  angle  of  attack 
of  a  wing  varies  with  time  as 
shown  in  sketch  (a),  the  resulting 
lift  coefficient,  also  shown  in 
sketch  (a),  is  designated  as  the 
indicial  lift  coefficient  due  to 
angle  of  attack.  Four  such  indi— 
cial  functions  will  be  evaluated. 


the  indicial  section  lift  coefficient  due  to  angle-of—  attack 
change  (without  pitching) 


the  indicial  section  pitching-moment  coefficient  due  to  angle— of— 
attack  change  (without  pitching)  measured  about  the  leading 
edge  and  considered  positive  when  the  trailing  edge  is  forced 
down 


the  indicial  section  lift  coefficient  due  to  pitching  on  a  wing 
rotating  about  its  leading  edge 


c  1  the  indicial  section  pitching-moment  coefficient  due  to  pitching 
^  measured  about  the  leading  edge  for  a  wing  rotating  about  its 

leading  edge  and  considered  positive  when  the  trailing  edge  is 
forced  down 


The  equations  which  transform  these  functions  to  those  for  a  wing 
pitching  about  a  point  a  distance  ac  back  from  the  leading  edge  and 
having  its  moment  center  a  distance  be  back  from  the  leading  edge  are 
simply 
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over  a  certain  planar  area  in  xzt 
space  representing  the  area  occupied 
by  the  wing,  and  elsewhere  pressure 
is  continuous.  For  a  wing  of  chord 
length  c,  this  area,  shown  shaded 
in  sketch  (h),  is  hounded  by  the 
traces  of  the  leading  and  trailing 
edges,  x  =  — V0tf  and  x  =  c— Vq^, 
and  the  line  along  which  the  motion 
started  ( i .  e . ,  the  x  axis ) . 

In  the  case  of  the  indicial 
function  due  to  pitch  for  a  wing 
rotating  about  its  leading  edge,  the 
boundary  condition  requires  that  the 
upwash  be  given  by  the  expression 

wu  =  -(x  +  yot«)  0  (6) 


(b) 


over  the  same  region  in  the  xt  plane  as  for  the  angle  of  attack  case 
and,  again,  that  pressure  be  continuous  elsewhere.  The  angle  of  pitch,  6 
is  taken  as  positive  when  the  trailing  edge  is  lower  than  the  leading 
edge,  and  §  is  the  time  derivative,  d0/dtf,  positive  when  the  trailing 
edge  is  falling  with  reference  to  the  leading  edge. 
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Direction  of  wing  motion 


Flight  path  of  leading 


(c) 


9  --  9  --  0 


The  difference  "between  8  and  a 
is  illustrated  in  sketch  (c).  The  angle 
of  attack  a  is  the  angle  between  the 
flat  wing  surface  and  the  tangent  to  the 
flight  path  of  the  leading  edge.  The 
angle  6  is  the  angle  between  the  flat 
wing  surface  and  the  horizontal. 

Sketch  (d)  shows  a  wing  undergoing  a 
sinusoidal  angle-of-attack  -variation 
with  a  zero  angle  of  pitch  throughout. 

For  convenience  and  in  order  to  distin¬ 
guish  from  the  pitching  wing,  such  a 
wing  will  be  referred  to  as  a  sinking 
wing.  Sketch  (e)  shows  a  wing  undergo¬ 
ing  a  sinusoidal  angle-of -pitch  varia¬ 
tion  taken  about  the  leading  edge  at  a 
constant  (zero)  angle  of  attack. 

The  variable  q  to  which  the  lift 
and  pitching-moment  coefficients  of  the 
pitching  wing  are  referred  is  equal  to 
0c /vQ,  a  dimensionless  form  of  expression 

for  rate  of  pitch. 


Adaptation  of  the  Indicial 
Functions  to  Maneuvers 


(d) 


When  the  wing  is  undergoing  a  maneu¬ 
ver  in  which  a  and  9  make  small  but 
arbitrary  deviations  from  zero,  the  bound¬ 
ary  condition  in  the  region  of  the  xt 
plane  traversed  by  the  wing  becomes 


wu  =  -VQ  a  -  (x  +  Y0t»)  e  (7) 


Since  the  theory  is  linear,  the  lift  and  moment  on  the  wing  are  calcu¬ 
lated  from  the  indicial  functions  by  the  principle  of  superposition. 
For  example,  the  lift  and  pitching-moment  coefficients  developed  by  an 
arbitrary  variation  of  a  and  8  are 


d 

Cl  dt7 


t» 


t* 


m 


dt’ 


c,  (t'-tx*)  a(tx* Jdti’  +  c,  *  (t’-tx’)  q(t1‘)dt1* 
_  6a 


(t’-tx*)  a(tx’)dtx‘  +  Cmq’  (t’-tx’ )  q(tx’)dtx’ 


>(8) 


o 


MCA  TN  2403 


9 


Equation  (8)  is  derived  for  an  angle— of-attack  variation  measured 
with  respect  to  the  flight  path  of  the  leading  edge  and  the  primes  on  the 
derivatives  indicate  that  the  wing  is  pitching  about  and  the  moments  are 
measured  with  respect  to  the  wing  leading  edge.  If  the  angle  of  attack 
is  measured  with  respect  to  the  flight  path  of  some  other  point  on  the 
wing  (as,  for  example,  the  center  of  gravity),  then  it  can  be  shown  for 
small  deflections  that 


cl  ~clx  — 


c£0 


<r 


where  ax  and  0i  are  the  angles  measured  with  respect  to  the  new 
flight  path  and  c£  is  the  distance  from  the  leading  edge  to  the  new 
reference  point.  (See  sketch  (c).)  The  total  value  of  the  section  lift 
and  pitching-moment  coefficients  for  an  arbitrary  variation  of  ai 
and  0!  would  be 


-akT 


al(ti* )  -  lqi(t!*) 


+  C2q,(-t,-ti,)qi(ta,)J-  dti* 


t» 


m 


’’H v/  {v(t,-tl’) 


a-iCti')  -  iqi(ti* )  j+  c  ,(t*-t1,)qi(ti,)T  liti* 

J  “q  J 


Adaptation  of  the  Indicial  Functions  to 
Flutter  Derivatives 


The  notation  adopted  in  this  report  does  not  coincide  with  the 
usual  notation  used  in  the  study  of  fluttering  wings;  however,  it  should 
be  sufficient  to  relate  the  present  notation  to  the  one  used  in  refer¬ 
ence  4.  Thus,  the  values  of  the  terms  Lx,  L2,  L3‘,  L4*  and  Mi’,  M2' , 
M3’ ,  M4‘  used  in  reference  4  can  be  evaluated  as  follows: 


Ia+iL2  =  -  _SL_ 

4k  t'-xadt’ 

.M1,+iM2*  = 


-iwt 


lim 


7  (t *— tx* )e1Utl*  dti‘ 
a 


.  —  iGdt 1 

ie _  lim  d 

2k  t >»  dt  * 


t» 


:>  (9) 


CjjjV  (t*— ti*  )eX  dtx* 


o 
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Lj+iL  *  = 


e— iwt'  4 


3  4  4k  2  t'^»4t 


rf 


c Za(t’—ti •  )+iu 


iuti' 

e  dti’ 


)(9) 


t* 


Cj^Ut’-tiM+i10  cm_  *  (4*— ti*  ) 


“q 


THE  SOLUTION  FOR  THE  INDICIAL  FUNCTIONS 


The  boundary  conditions  which  represent  a  wing  undergoing  a  sudden 
jump  in  pitching  velocity  or  angle  of  attack  have  been  discussed  and 
presented  in  equations  (5)  and  (6).  Their  application  to  wings  flying 
at  several  Mach  numbers  will  now  be  presented.  The  incompressible  case 
will  first  be  discussed,  and  then  the  cases  in  which  Mq  equals  0.8, 
1.0,  1.2,  and  2.0  will  be  considered. 


Incompressible  Case,  MQ  =  0 


The  solution  for  the  incompressible  case  is  applicable  when  the 
forward  speed  of  the  wing  is  small  compared  with  the  speed  of  sound  so 
that  the  ratio  V0/a0  can  be  neglected  in  comparison  to  unity.  The 
basic  partial  differential  equation  governing  the  flow  field  was  presented 
in  equation  (2)  as 

^xx  +  ^zz  =  0 

subject  to  the  boundary  conditions  discussed,  and  the  equation  for  the 
loading  coefficient  was  given  in  equation  (4)  as 

=  _k_ 

°  V02  dt' 

This  boundary— value  problem  corresponds  precisely  to  that  which  is 
studied  in  three-dimensional,  steady-state,  lifting-surface  theory  under 
the  classification  "s  lender-wing  theory."  (See  references  6  and  7.) 

Such  an  analogy  is  useful  since  well-established  concepts  in  one  theory 
can  be  immediately  carried  over  into  the  other,  (it  should  be  re¬ 
emphasized,  perhaps,  that  the  subsequent  treatment  of  the  incompressible 
case  is  not  intended  to  be  an  improvement  on  Wagner fs  original  deriva¬ 
tion  but  rather  it  is  a  rederivation  along  lines  that  will  be  used  later 
in  the  analysis  of  the  compressible  cases.) 
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The  initial  pulse .—  The  first  analogy 
with  slender-wing  theory  which  will  he  used 
concerns  the  initial  pulses  that  occur  in  the 
values  of  lift  and  pitching  moment.  It  is  a 
well-known  result  (reference  6)  that  the 
total  lift,  as  given  hy  slender-wing  theory, 
on  the  wing  shown  in  sketch  (f)  is  a  func¬ 
tion  only  of  the  maximum  span  and  the  value 
of  wu  along  the  section  of  maximum  span 
(section  AA).  It  is,  therefore,  independent 
of  the  wing  twist  and  leading-edge  shape 
ahead  of  section  AA.  This  concept  has  been 
extended  in  slender-wing  theory  to  the  extreme 
case  shown  in  sketch  (g)  of  a  rectangular 
wing.  The  lift  on  such  a  wing  is  concen¬ 
trated  entirely  along  the  leading  edge  and  is 
a  function  only  of  the  span  of  that  edge  and 
the  value  of  wu  there.  By  the  analogy 
existing  between  the  two  theories,  therefore, 
it  is  evident  that  the  solution  to  the  indi— 
cial  problems  in  two-dimensional,  incompress¬ 
ible,  unsteady  flow  (sketch  (h))  will  contain 
a  pulse  at  t*=0. 

The  evaluation  of  this  pulse  will  be 
treated  briefly.  A  solution  to  equa¬ 
tion  (2)  for  the  vertical  induced  velocity 
in  the  z=0  plane  can  be  written  in  terms 
of  the  jump  in  u  across  the  z=0  plane 
(see  reference  7)>  thus,  for  the  shaded 
area  in  sketch  (i)  this  is 

“»<*>  -  -  Sf  f  (10) 

u-a. 

The  general  inversion  of  equation  (10) 
can  he  written 


Pulse  in 
lift  along 
leading 
edge 

(9) 


Pulse  in  lift  at  t'=  O. 
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In  the  present  case  A  is  zero,  since  A<p  is  zero  at  x=— a  and  x=b, 
and  an  integration  of  both  sides  of  equation  (11)  with  respect  to  x 
between  the  limits  b  and  x  gives 


Acp(t',x) 


w(t»,Xi) 


In  "L  Cb-SiHa+x)  +  J (a+Xi )  (b— x ) 

x-xi)(a+b) 


Adoption  of  the  notation 


v(t,,x1)  +  constant  = 


w(t,,x1)  dxi 


and  integration  by  parts  leads  to  the  equation2 

Acp  -  _  2«/(b-x)(a+x)  p  V(t«,xi)  dxj.  . 

rt  J-a  (MiV(b-Xi )  (a+Xi ) 

The  loading  can  now  be  determined  by  using  equation  (4).  If  the 
shaded  area  in  sketch  (i)  is  allowed  to  vanish,  all  the  loading  accumu¬ 
lates  along  the  x  axis  in  the  region  0<  x<  c.  Therefore,  the  inte¬ 
gral  of  the  loading  with  respect  to-  t‘  over  the  shaded  portion  must  be 
considered.  The  final  result  for  the  pulse  loading  (Ap/qQ)  at  t*=0 

can  be  expressed  in  terms  of  the  6  function  (see  list  of  symbols)  as 


rC  v(0,Xi)  dxx 

J  (x-Xi  )  </(  C-Xj.)  Xi 
o 


(13) 


The  boundary  conditions  for  the  sinking  and  pitching  wing  given  by 
equations  (5)  and  (6),  when  inserted  into  equation  (13),  yield 


The  constant  gives  zero  when  placed  in  equation  (12),  provided 
h  >  x  >  -  a. 
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After  integration,  the  pulse  -values  for  lift  and  pitching  moment 
may  also  he  obtained.  Hence, 


*r6(0) 


SO  --sC6(0) 

o 

where  the  primes  indicate  that  the  wings  are  pitching  about  and  the 
moments  are  measured  about  the  leading  edge.  These  expressions  may  be 
inserted  in  equations  (8)  and  (9)  and,  since  the  integrand  becomes  zero 
everywhere  except  at  the  point  tif=tf,  the  expressions  for  the  lift  and 
pitching-moment  coefficient  developed  by.  an  arbitrary  variation  of  a 
and  0  with  time  are 


a  +  6 
2Y  kY  2 

o  o 


9_  jtcf  'Q 
64  v  * 


The  variation  for  t*  >  0.-  The  integral  equation  (10)  is  still  per¬ 
fectly  valid  when  applied  to  the  flow  field  for  t*  >  0.  It  is  convenient 
to  rewrite  the  equation  in  this  case,  however,  so  that  the  effects  of  the 
vorticity  on  the  wing  and  in  the  wake  are  separated.  Thus, 


wu  = 


2*  J 


’c-Srot’  Au(t«,Xl) 

-  dxi  — 

X-X;l 


j 

-v  t* 

o 


Au* (xi ) 


tei  (l6) 


C-V’ 
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where  Au*(xx)  is  the  value  of  Au  in  the 
starting  vortex  wake.  It  is  independent 
of  t'  since  its  value  at  all  points  along 
the  line  ab  in  sketch  (j)  is  the  same  as 
at  the  point  a. 

A  reduction  of  equation  (l6)  can  be 
obtained  for  the  case  of  the  sinking  wing, 
where  wu=  — YQa,  by  using  the  inversion 
given  by  equation  (ll).  Thus, 


JtAu(t*,x)  J ( x+Y Qt * ) ( c-V Qt *-x )  =  A-tf(2x+2Y  t’-c)  Ya 


r°  Ti  * 

J  L  (x-Xx) 


c— V  Qt  * 


dx 


Since  A  is  given  by  the  relation 


A  =  /  AuCt’^x)  dxx  =  -  / 

-V  c— V  Qt  * 


Au*(xx)  dxx 


it  follows  that 


.  , 


7(x+Y0t»)(c-V0t'-x)  «  y(x+V0t»)(c-Vot'-x) 


c-Y  t» 
o 


—  j (xi+V0t  * )  (xx-c+Y0t  • )  dxx 


(17) 


According  to  the  Kutta  condition  Au(t',x)  vanishes  as  x  approaches 
C-V>  integral  equation  for  Au*(xx)  thus  becomes 


Voa  =  - 


Au* (xx ) 


r 

h 


Xl+Vpt  * 
Xx-c+V0t 1 


dxx 


(18) 


c— Y Qt  * 
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which  was  derived  and  studied  originally  toy  Wagner  (reference  1). 

The  section  lift  and  pitching  moment  can  toe  derived  in  terms  of 
Au* ( x x )  in  the  following  manner.  By  definition,  the  section  lift  l 
is 


l 


C— \/A  t  ’ 


PnV, 


O’O 


-V’ 


! 


c-Y0t* 

°  &CD 


-V’ 


dx 


(19) 


Since  the  value  of  Acp  is  zero  at  the  leading  edge  and  at  the  trailing 
edge  is  equal  to  the  total  circulation  T,  two  alternative  forms  for  the 
lift  can  he  written 


a  rtf 

1  =  PoV  +  po~^T J  ^p(t*,x)  dx 


and 


-V’ 


c-V0t* 


1  '  po(c-V)  o  to 


(20) 


(21) 


By  substituting  equation  (17)  into  (21)  and  integrating,  it  can  toe  shown 
that 


l  = 


T r  2  PoCVo 
rtp0cV0  a  +  ° 


c-V0t» 


Au*(xi) 

'  —  -  - =  clXi 

V  (xi+Y0t,)(x1-c+Y0t’) 


(22) 


Since  the  value  of  Au*(xx)  has 
been  determined  toy  Wagner,  equa¬ 
tion  (22)  represents  a  solution 
for  the  section  lift.  A  plot 
of  its  variation  in  coefficient 
form  is  shown  in  sketch  (k). 
Initially  there  is  the  pulse 
having  an  intensity  defined  toy 
equation  (15).  After  the  pulse 
at  t*=0,  the  value  of  the  sec¬ 
tion  lift  coefficient  starts  at 
one-half  its  asymptotic  value. 

It  then  increases,  slowly 
approaching  its  asymptote  of  2rt. 


(k) 


V0t‘/c 
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By  definition,  the  pitching  moment  can  be  written 


c-V0t* 

m  =  -5poV/  (T0t'«)&ai 

-V  t' 


rc-V' 

'  -  p°77  t,  (V*«>  ^ 


(23) 


where  the  moment  is  taken  about  the  leading  edge  and  is  positive  when 
the  trailing  edge  is  forced  down.  A  development,  similar  to  the  one 
given  for  the  lift,  gives 


m 


(24) 


This  result,  that  for  t  >0  the  indicial  center  of  pressure  remains  con¬ 
stant  at  the  quarter  chord  throughout  the  motion,  is  classical. 

If  the  boundary  condition  for  a  pitching  wing,  wu  =  -(x+V0t’)  0, 
is  substituted  into  equation  (l6)  and  the  inversion  given  by  equation  (ll) 
is  again  used,  it  can  be  shown  in  the  same  manner  used  in  the  derivation 
of  equation  (l8)  that  for  x  =  c— V0t’  the  relation 


3 

¥ 


Au*(xi) 


Xi+Vpt' 
Xi~c+V  0t f 


dxi 


applies.  This  integral  equation  applies  to  a  wing  pitching  about  its 
leading  edge.  If,  instead,  the  wing  is  pitching  about  the  three- 
quarter-chord  position,  an  essential  simplification  is  achieved.  In  this 
latter  case,  downwash  is  given  by  the  expression 


wu  =  -0(x+Vot‘-|  c) 


(25) 


and  the  resulting  integral  equation  becomes 

°  =  r  £u  „*(xi) 


h+U' 


c-V 


Xi— c+V_t 


T 


(26) 
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% 


A 


where  Au3*(xi)  represents  the  vorticity  in  the  wake  following  such  a 
motion.  The  solution  to  equation  (26)  is  simply 

Au3*  (xi)  =  0  (27) 


From  equation  (27)  it  follows  that  the  total  indicial  lift  for  t*>0 
on  a  wing  pitching  about  the  three-quarter  chord  point  is  zero,  and  that 
the  wing  wake  is  free  of  vorticity.  Further,  it  can  be  shown  that  the 
total  indicial  pitching  moment  (still  measured  about  the  leading  edge) 
is 


pitc*3 

m  =  "^ 


V0  & 


The  transfer  of  equations  (27)  and  (28)  hack  to  the  case  in  which 
the  wing  is  pitching  about  its  leading  edge  can  be  readily  accomplished 
by  means  of  the  boundary  conditions  shown  in  sketch  (l).  Hence,  if 
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By  means  of  equations  (24) 
indicial  functions,  cma’ ,  ci^J , 
the  indicial  lift  function  for 


and  (29)  the  expressions  for  the  three 
and  Cjn^’  can  all  be  written  in  terms  of 
t  >  0 .  Hence , 


Cm 


a 


c 


c 


l 


a 


cV 


_3  „  _£ 

16  8 


>  (30) 


The  variations  of  the  four  indicial  functions  are  shown  in  figure  1. 
For  values  of  Y0t*/c  larger  than  those  shown  in  the  figure  the  approxi¬ 
mate  equation  suggested  in  reference  2  can  be  used,  namely. 


c 


a 


2it 


This  alternative  result  has,  according  to  reference  2,  an  error  of  2  per¬ 
cent  or  less  for  the  entire  range  of  time  from  0+  to  infinity. 


Subsonic  Case,  MQ  =  0.8 


When  the  Mach  number  is  no  longer  small,  the  analysis  in  the  pre¬ 
ceding  section  must  be  modified.  As  an  example  of  this  modification,  the 
calculations  for  a  wing  traveling  at  a  Mach  number  equal  to  0.8  will  be 
carried  out  in  detail. 

Equation  (3)  presents  the  basic  partial  differential  equation  of 
the  flow  field,  and  equation  (^)  gives  the  expression  for  the  loading 
coefficient.  The  analogy  which  existed  in  the  incompressible  case 
between  the  theory  for  the  unsteady,  two-dimensional  wing  and  slender— 
wing  theory  exists  in  this  case  between  the  theory  for  the  unsteady, 
two-dimensional  wing  and  the  theory  for  a  steady— state,  three-dimensional 
wing  traveling  at  a  supersonic  speed.  Thus,  in  the  three-dimensional, 
steady-state  case  the  partial  differential  equation  governing  the  flow  is 

<Pyy  +  ?zz  =  P2(Pxx  (31) 

and  the  equation  for  the  loading  coefficient  is 
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AP_.j-.dcp 
<1  Y0  Sx 


(32) 


The  boundary  conditions  are  in  both  cases  that  epz  is  given  over  a 
portion  of  the  z=0  plane.  It  is  evident  by  a  comparison  of  equa¬ 
tions  (3)  and  (31)  and  equations  (4)  and  (32)  that  results  from  the 
three-dimensional,  supersonic,  steady— state  study  (hereinafter  referred 
to  as  the  steady-state  case)  can  be  transferred  to  the  two-dimensional, 
unsteady  study  (hereinafter  referred  to  as  the  unsteady  case)  simply  by 
replacing  x,  y,  and  3  in  the  former  case  by  t,  x,  and  1,  respectively, 
and  by  dividing  the  result  for  the 
loading  coefficient  by  Mq. 


The  analog  to  the  boundary 
condition  for  the  problem  of 
finding  the  indicial  loading  on  a 
two  dimensional  wing  flying  at  a 
Mach  number  equal  to  0.8 
(sketch  (m))  is  the  boundary  con¬ 
dition  for  the  problem  of  finding 
the  loading  on  a  constant-chord, 
swept-forward  wing  tip  with  a 
subsonic  trailing  edge  such  as 
that  shown  in  sketch  (n).  The 
Mach  cones  in  the  steady-state 
case,  traces  of  which  are  shown 
as  dotted  lines  in  sketch  (n), 
become,  in  the  unsteady— state 
analog,  the  locus  of  the  sound 
waves  which  started  at  t=0  from 
the  leading  and  trailing  edges  of 
the  two-dimensional  wing 
(sketch  (m)).  Finally,  the 
analog  in  the  steady— state  field 
of  the  unsteady  wing  would  be  a 
flat  plate  for  the  unsteady 
sinking  wing  and  a  plate  with  a 
linear  variation  of  twist  for 
the  unsteady  pitching  wing. 

A  detailed  analysis  of  an 
unsteady,  two-dimensional  wing 
flying  at  a  Mach  number  equal  to 
0.8  will  be  presented.  Just  as 
in  the  section  on  incompressible 
flew,  the  study  will  be  divided 
into  two  parts.  In  cases  for 
which  Mq^O,  however,  the  indi— 
cial  functions  contain  no  pulse 
at  t=0.  Hence,  the  first  part 
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of  the  Mo=0.8  study  will  be  concerned  with  the  behavior  of  the  indicial 
functions  in  an  interval  for  which  t  is  small  but  finite  (actually 
0  <  t  <  5c)  and  the  second  part,  with  their  asymptotic  behavior. 

The  early  stage.-  The  analog  which  exists  between  the  steady-state 
and  unsteady  cases  may  be  utilized  to  great  advantage  since  the  large 
number  of  special  methods  and  techniques  developed  for  the  solution  of 
problems  in  the  former  case  may  be  applied  to  the  solution  of  the 
analogous  problems  in  the  latter  field.  In  this  manner  an  exact  solu¬ 
tion  for  the  loading  over  the  five  regions  shown  in  sketch  (o)  may  be 
obtained  for  both  the  sinking  and  the  pitching  wing  in  an  unsteady  flow 

field  by  the  use  of  methods  such 


as  those  presented  in  reference  8. 
Solutions  for  larger  values  of 
time  could  also  be  obtained,  but 
the  labor  involved  in  calculating 
such  cases  becomes  prohibitive  and, 
as  will  be  shown  later,  approx¬ 
imate  methods  can  be  developed 
which  extend  the  solutions  for  the 
indicial  lift  and  pitching— moment 
curves  to  their  asymptotic  values. 


X  -  “  M0ty/ 

'-V7 


% 


!  +  Mo 


/  *  ^ 
V  / 


I -Mo 


The  analysis  used  to  calculate 

(o)  the  loading  in  terms  of  xD  =  x/c 

and  to  =  t/c  over  the  five 

regions  shown  in  sketch  (o)  on  both  the  sinking  and  pitching  wings  is 
outlined  in  appendix  A.  Plots  of  the  loading  on  a  sinking  wing  are 
shown  in  sketch  (p)  and,  in  more  detail,  together  with  the  loading  on  a 

pitching  wing,  in  figure  2.  At  to=0  the 
api  >  loading  is  constant,3  and  as  time  increases 

the  loading— coefficient  curve  approaches  the 
—  familiar  two-dimensional,  steady— state  shape 

«  tV^ST"  given,  for  the  sinking  wing,  by  the  equation 


Ap  _  4a  I  c— x 
qQ  3  v  x 


The  result  that  the  initial  shape  of  the  load  distribution  is  the  same 
as  the  shape  of  the  given  curve  for  wu  also  applies  to  all  three- 
dimensional  wings  of  arbitrary  plan  forms  traveling  at  subsonic  or 
supersonic  Mach  numbers . 
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where  x  is  the  distance  from  the  leading  edge  which  is  at  the  point 
x=0.  Sketch  (q)  shows  that  the  distribution  at  t0  =  t/c  =  l/(l-Mo)  =  5 
is  already  essentially  the  same  distribution  as  that  obtained  at  t0=co 
(i.e.,  the  agreement  is  good  with  the  curve  produced  by  multiplying  the 
right  side  of  equation  (33)  by  a  constant  factor).4  The  use  of  this  fact 
simplifies  subsequent  analysis  concerning  the  asymptotic  behavior  of  the 
indicial  lift  curve. 


The  indicial  lift  and  pitching-moment  functions  were  also  calculated 
(see  appendix  A)  in  the  range  0  <  to  l/(l-Mo).  Their  variation  in 
this  interval  is  shown  in  sketch  (r)  for  Mq  =  0.8.  It  is  evident  that 
the  calculations  must  be  extended  beyond  the  point  t0=5  since  the 
asymptotic  values  are  not  yet  even  closely  approached. 

Before  studying  the  nature  of  these  curves  for  large  values  of  tQ, 
however,  it  is  useful  to  examine  them  with  reference  to  the  discussion 
in  the  previous  section  on  incompressible  flow.  Bor  example,  it  was 
pointed  out  that  the  indicial  center  of  pressure  on  the  sinking  wing 
remained  at  the  quarter-chord  point  for  tT>0.  It  is,  therefore,  per¬ 
tinent  to  consider  the  location  of  the  center  of  pressure  on  the  sinking 
wing  when  Mq  =  0.8.  By  means  of  the  curves  given  for  c}  and  c™  f 
and  by  the  relationship  a  ^ 


A  similar  result  was  noted  in  the  study  of  the  load  distribution  on 
swept— back  wings  with  subsonic  leading  edges  (reference  9)* 
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the  variation  of  (x/c)c^  is  easily  evaluated  ((x/c)c#p#  is  the  dis¬ 
tance  between  t'he  leading  edge  and  the 'center  of  pressure  divided  by  the 
total  wing  chord).  This  variation  is  shown  in  sketch  (s).  It  is 
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apparent  that  the  center  of  pressure  is  very  close  to  the  quarter-chord 
point  for  values  of  t0>5.  Hence,  at  Mq  =  0.8  the  significant  effect 
of  compressibility  on  the  location  of  the  center  of  pressure  is  contained 
in  the  interval  0  <  to  <  5.  This  leads  immediately  to  the  further  use¬ 
ful  conclusion  that  for  values  of  tQ  greater  than  5  the  value  of  the 
pitching— moment  coefficient  is  given  essentially  by  the  equation 

Cma*  =  "  cZa/4 

Likewise,  it  is  apparent  from  the  discussion  of  the  incompressible 
case  that  the  indicial  functions  for  the  pitching  wing  can  also  be 


0 
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expressed  in  a  more  convenient  form  by  shifting  the  axis  of  rotation 
from  the  leading  edge  to  the  three-quarter-chord  point.  Using  the  values 
of  c 2a,  cn\x* >  clq*  >  an^  Cmq*  Presented  in  sketch  (r),  the  variation  of 

(c^)3/4  and  (cm^f  )0/-  may  be  calculated  from  the  definitions  given  in 


(Cmq*  ^3/4 


equation  (29).  These  curves  are  presented  in  sketches  (t)  and  (u)  and 
again  it  is  apparent  that  at  a  Mach  number  equal  to  0.8  the  compressi¬ 
bility  effects  are  contained  in  the  interval  0  <  t0  <  5. 

The  later  stage.—  It  follows  from  the  preceding  discussion  that  when 
tQ  is  large,  the  values  of  the  indicial  functions  cm* ,  c7  * ,  and  cm  * 

for  compressible  flow  can  also  be  expressed  adequately  in  terms  of  c7a 
by  equations  similar  to  equations  (30)  which  were  derived  for  incompress¬ 
ible  flow.  Thus  for  tQ  >  5,  one  can  write 

cma’  =  “  cZaA 


la'  =  3 


cn^*  =  “  (3cla/l6)  -  (k/80) 

The  four  indicial  responses  have  all  been  shown  to  depend  only  on 

the  value  of  the  function  c-i  if  tQ  is  large.  It  remains,  therefore, 

La 

to  determine  the  asymptotic  behavior  of  c^  •  This  can  be  accomplished 
in  the  following  way.  Consider  the  steady-state  solution  for  the  lift  on 
a  two-dimensional,  flat  lifting  surface  traveling  at  a  subsonic  Mach 
number.  As  was  pointed  out  by  Wieghardt  (reference  10),  if  the  lift  on 
such  a  surface  is  represented  by  placing  at  the  quarter— chord  point  a 
vortex  which  has  the  same  circulation  as  that  developed  by  the  wing,  the 
angle  of  attack  measured  at  the  three— quarter— chord  point  will  be  the 
same  as  that  of  the  simulated  lifting  surface.  Extending  this  concept  to 
include  the  unsteady  effects,  an  investigation  will  be  made  of  the  varia¬ 
tion  with  time  of  the  vortex  strength  which  will  maintain  a  constant 
angle  of  attack  at  the  three— quarter— chord  station  following  an  impulsive 
start  at  to=0. 

The  analogous  problem  in  steady— state  theory  becomes  one  of  finding 
the  strength  of  the  vortex  system,  shown  in  sketch  (v),  which  gives  a 
constant  value  of  w  along  the  line  CD.5 

Each  vortex  compos ing  this  system  lies  along  the  line  AB,  extending 
from  the  minus  infinity  toward  the  origin,  and  trails  back  parallel  to 


In  the  vicinity  of  the  origin,  of  course,  this  representation  gives  a 
poor  approximation  to  the  original  boundary— value  problem;  hence,  use 
of  the  results  must  be  limited  to  the  regions  in  which  to  is  large. 
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(v)  Steady  -  state  analog  Unsteady  case 


the  Xi  axis  to  form  the  trailing  vortex  sheet.  Note  that,  for  con¬ 
venience,  the  origin  of  the  axis  system  has  been  located  at  the- quarter- 
chord  point.  The  solution  to  such  a  problem  in  steady— state,  lifting- 
line  theory  would  result  from  the  solution  of  the  integral  equation6 
(for  a  development  using  the  notation  adopted  here,  see  reference  11): 


_ A<P(yi)  dxi _ 

[(«x)a  -  02(y-yi)233/2 


where  A9  is  not  a  function  of  Xi  since  the  strength  of  a  trailing 
vortex  is,  of  course,  constant.  The  area  of  integration  t  is  the 
region  within  the  forecone  springing  from  the  point  x,y. 

If  the  above  equation  is  transformed  by  means  of  the  analogy  to 
represent  the  solution  of  the  unsteady  problem  (see  sketch  (w)), 

e 

The  symbols 

be  taken.  Thus  (see  reference  11  or  12), 


are  used  to  indicate  that  the  finite  part  is  to 


dx 


(x,b)  -  G 


(x,a) 


where  G  (x,y)  is  the  indefinite  integral  of  f(y)/(x— y)2.  Further, 


(x-y)3/2 


»-e 

5x 


dy 
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3,  x,  and  y  are  replaced  by  1,  t,  and  x,  respectively;  and  Aq>,  the 
total  jump  in  potential  at  a  given  section,  is  replaced  by  the  circula¬ 
tion  T.  Hence. 

M  =  jf  tol  /  dti 

where  T,  as  indicated  in  sketch  (w),  is  the  area  in  the  forecone  from 
the  point  P  which  lies  always  along  the  linfe  x  =  (l/2)c  -  Mot. 
Integration  with  respect  to  ti  reduces  the  last  equation  to 


w 


i/° 

J  Mq(x— t) 


1  +  Mq 


which,  by  means  of  the  substitution  xx/c  =  x2/2  becomes  along  the  line 
x  =  |  -  Mot 


1  ( 40+n q— x2 )  r(x2) 

0Jt3  X  (n2-X0+x2)2  J (X0-x2)(Xq-Hi1-x2) 


(35) 


where  Xq  -  2  Moto  -  do*  do  -  Mo/(l+Mo),  Hi  -  2  Mq/(1— Mq2),  |i2  =  l/(l+Mo), 
tQ  =  t/c,  and  where,  of  course,  w  is  a  constant  equal  to  — VQa. 

A  solution  for  r(x2)  in  the  integral  equation  (35)  may  be  obtained 
by  expanding  r  in  a  series  of  the  form 


r(x2) 


rtcaV0 

“T““ 


X2 

( a0+x2 ) 


+  bi 


xg 

(a0+x2)3 


+ 


(36) 


Place  equation  (36)  into  (35)  and  expand  in  powers  of  l/X0.  There 
results  the  expression 

_ W_  =  .  £1  Cg2n(l/x0)  _C3_ 

VQa  X  X0  \02  X03  +  '  *  * 


in  which 


MCA  TN  2403 


27 


Ci  -  t>i 


.  ao  .  1 
2  02 


Hence,  if  aQ  and  bi  are  chosen  so  that  Ci  is  equal  to  zero,  an 
expression  for  T  will  he  obtained  which  represents  the  solution  to  the 
integral  equation  (35)  correct  to  the  first  order  in  l/tQ  (i.e.,  l/XQ) 
for  large  values  of  t0.  Further,  if  equation  (36)  is  expanded  in 
powers  of  l/x2,  there  results 


1  +  — 7p)  +  •  •  • 

which  becomes,  using  the  condition  for  Ci  and  relating  x2  and  tQ 
by  the  equation  of  the  leading  edge. 


rfe)  =  IS*. 


r(t0) 


jtca,V0 

~T~ 


2Mot0P2 


(37) 


The  relation  between  circulation  and  lift  in  incompressible  flow 
has  been  derived  and  presented  as  equation  (20).  For  compressible  flow, 
this  expression  becomes 


n  nl— M oto 

1  =  P0V0r  +  Poao  T7~  /  A<p(to,x0)dx0 

6t°  -'Moto 

where  xQ  =  x/c.  In  order  to  obtain  a  complete  expression  for  the  sec¬ 
tion  lift,  it  is  obviously  necessary  to  know  the  chordwise  variation  of 
A9.  The  asymptotic  behavior  of  l  can  be  calculated,  however,  by 
applying  the  results  presented  in  sketch  (q).  This  result  suggests  that 
for  large  values  of  tQ  the  value  of  A<P(t0,x0)  used  in  the  equation 
for  section  lift  can  be  expressed  by  the  product  of  Ap(oo,x0)  anil 
r('to)/r(“)«  In  other  words,  for  large  values  of  tQ  the  shape  of  the 
chordwise  distribution  of  vorticity  is  the  same  as  the  two-dimensional, 
steady— state  value.  An  indication  of  the  accuracy  of  such  an  approxima¬ 
tion  is  shown  in  sketch  (x)  where  the  precise  value  of  A*?(t0,x0)  is 
compared  with  the  approximation  at  tc=5. 

Since  T(oo)  =  jtaV0c/3,  the  substitution  of 

AtKt0,x0)  =  Ap ( 00, xq )  r(t0)/r(<») 


1  =  p0vQr  +  l  Poao 


gives  for  the  lift 
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Place  equation  (37)  in  the  above  equation  and  it  is  evident  that,  to  the 
first  order  in  l/t0,  the  section  lift  is  given  by  p0VoT.  Hence,  it 
follows  that  for  large  t© 

=  f  (  1  “  2M0t0P2  +  *  '  ')  (38a) 

There  remains  the  problem  of  joining  the  above  result  for  c ia  with 
the  one  derived  in  the  preceding  section  and  valid  for  0<to<5*  To 
accomplish  this  end  the  equation 


■a 


2* 

P 


1  - 


hi 


h0  +  2Mot0P2  (h0  +  2MotoP2)2 


(38b) 


was  used  to  express  c->  for  the  range  5^t0.  Obviously  the  value  of 
cZ(X  given  by  equation  (38b)  has,  to  the  first  order  in  l/t0,  the  same 

asymptotic  variation  as  that  given  by  equation  (38a)  regardless  of  the 
values  of  the  constants  hD  and  hi.  These  constants  can  be  chosen, 
therefore,  so  that  both  the  magnitude  and  slope  of  the  indicial  ci 
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curve  given  by  equation  (38b)  are  continuous  at  t0=5  with  the  exact 
c^a  curve  obtained  for  0  <  to  ^  5  in  the  previous  section. 


The  curves  for  c^*,  c^1 ,  and  c^1,  as  previously  discussed,  can 
be  calculated  from  cj^  by  the  use  of  equations  (34).  In  each  of  these 
three  cases,  the  magnitude  of  hi  was  modified  slightly  so  the  resulting 
curve  would  be  continuous  at  t0=5  with  the  exact  results  presented  in 
the  preceding  section.  The  final  expressions  are 


_  2jt 

C^a  3 


t 


— 7t 


%  ~  20 


C7  *  =  32. 

l<l  23 


1  - 


—it 


"  23 


1.736 

131.2  1 

17.06  +  to 

(17.06  +  to)2  J 

1.736 

121.8 

17.06  +  t0 

(17.O6  +  t0)2 

1.736 

134.3 

17.06  +  to 

(17.06  +  t0)2  . 

1.302 

90.53 

17.06  +  to 

(17.06  +  t0)2 

>  (39) 


J 


The  final  indie ial  section  lift  and  pitching— moment  curves  are 
shown  in  figure  1  plotted  against  the . parameter  V0tf/c,  the  number  of 
chord  lengths  traveled  (V0t*/c  =  Mot0).  Tabular  results  of  the  indi- 
cial  curves  are  also  presented  in  table  I.  ' 


Sonic  Case,  Mq  =  1.0 


The  general  results,  obtained  in  the  preceding  section  and  pre¬ 
sented  in  appendix  A,  for  the  indicial  loading  over  the  sinking  and 
pitching  wing  may  be  extended  to  the  sonic  case.  Furthermore,  the  two 
intervals  for  which  analytic  results  in  a  closed  form  were  presented  in 
appendix  A  now  cover  the  complete  time  range  since  0  <  tQ  <  l/(l+Mo) 
becomes  0  £  t0  ^  0*5  and  l/(l+Mb)  ^  tQ  <  l/(l-Mo)  becomes 
0.5  <  tQ  -  Hence,  by  an  appropriate  limiting  process  equations  (A8), 
(A9),  (A10),  and  (All)  become  for  0  ^  tQ  ^  0.5 

c^a  =  h 


Cmrx* - 2  +  tQ2 

c =  2  +  to2 

cmq’  =  -(4/3)  -  (2/3)t03 


(40a) 
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and  for  0.5  ^  tQ  ^  00 


Since  the  magnitude  of  the  functions  in  equation  (40b)  grows  indef¬ 
initely  with  increasing  time,  the  assumptions  of  linear  theory  are  even¬ 
tually  violated.  However,  for  moderate  values  of  tQ,  these  functions 
have  the  same  order  of  magnitude  as  similar  indie ial  curves  for  Mach 
numbers  other  than  1.  These  effects  are  illustrated  in  figure  1. 


* 


Supersonic  Case,  Mq  =  1.2  and  2 

The  method  of  obtaining  solutions  for  the  indicial  functions  at 
supersonic  Mach  numbers  parallels  the  development  presented  for  Mq  =  0.8. 
The  steady-state  analog  to  the  supersonic  unsteady  wing  problem  is  a 
constant-chord  wing  tip  with  a  supersonic  trailing  edge.  See  sketches  (y) 
and  (z).  It  is  well  known  that  the  problem  of  finding  the  loading  over 
wing  plan  forms  with  all  supersonic  edges  is  one  of  the  simplest  in  three- 
dimensional,  lifting— surface  theory.  In  fact,  since  the  upper  and  lower 
surfaces  are  noninteracting,  the  solution  is  determined  by  integrating 
sources  within  the  Mach  forecone.  The  analysis  for  C2a  has  already 
been  carried  out  in  reference  5* 
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The  analysis  used  to  calculate  the  loading  in  terms  of  x0=x/c  nnfl 
t0=t/c  over  the  three  regions  shown  in  sketch  (aa)  is  outlined  in 
appendix  B.  An  example  of  the  manner  in  which  the  loading  varies  with 
time  over  a  sinking  wing  traveling  at  a  Mach  number  equal  to  1.2  is  given 
in  sketch  (ab)  and,  in  more  detail,  together  with  the  loading  over  a 
pitching  wing,  in  figure  3. 


The  expressions  for  the  indicial  lift  and  pitching— moment  coeffi¬ 
cients  are  given  analytically  in  appendix  B,  and  plotted  in  figure  1. 

It  can  be  shown  that  the  results  given  in  appendix  B  reduce  to  the 
expressions  given  by  equations  (40)  when  Mq  is  allowed  to  approach  one, 
so  that  there  is  no  discontinuity  in  the  theory  in  passing  through  the 
sonic  range. 


CONCLUDING  REMARKS 


The  use  of  the  analogy  between  the  basic  flow  equations  in  steady- 
state,  lifting-surface  theory  and  in  unsteady-state,  airfoil  theory  has 
resulted  in  a  method  of  calculating  two-dimensional  indicial  functions 
throughout  the  subsonic  and  supersonic  flight  range.  The  results  are, 
of  course,  subject  to  the  restrictions  of  linearized,  compressible-flow 
theory  and,  for  example,  the  calculated  responses  given  in  figure  1  for 
sonic  speeds  must  be  considered  as  being  outside  the  realm  of  validity 
within  a  few  chord  lengths  of  travel.  In  application  to  high-frequency 
oscillations,  however,  the  initial  portions  of  the  indicial  curves 
dominate  the  response  characteristics  of  the  airfoil  and  calculations 
near  Mq  equal  to  one  need  not  be  invalid. 

In  the  supersonic  range,  the  expressions  in  appendix  B  are  calcu¬ 
lated  for  arbitrary  Mq.  It  is  apparent  that  no  fixed  Mach  number  effect 
can  be  used  in  transient  responses  except  at  high  values  of  flight  Mach 
number.  In  the  subsonic  range,  however,  the  expressions  apply  only  for 
the  period  of  time  0  £  t0  £  I/CI-Mq).  For  values  of  tQ  i  l/(l-Mo), 
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the  method  outlined  in  the  report  appears  to  be  satisfactory  for  Mach 
numbers  equal  to  or  greater  than  0,8.  Preliminary  calculations  made  at 
Mach  numbers  less  than  0.8  indicate  that  it  is  necessary  to  extend  the 
exact  solutions  past  t0  =  l/(l-Mo).  This  extension  is  feasible  if  the 
cancellation  techniques  outlined  in  reference  13  are  employed  and  the 
more  difficult  integrals  are  expanded  into  series  form. 

Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 

Moffett  Field,  Calif.,  April  12,  1951. 
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APPENDIX  A 

DETERMINATION  CF  SUBSONIC,  INDICIAL,  SECTION  LIFT  AND 
PITCHING-MOMENT  CURVES 
THE  LOAD  DISTRIBUTION 


The  following  results  for  the  indicial  load  distribution  on  sinking 
or  pitching  wings  can  be  obtained  in  two  ways.  One  of  these  methods  will 
be  outlined  in  the  subsequent  paragraphs.  The  other  is  outlined  in 
references  13  and  1^  and  is  referred  to  as  the  lift-cancellation  techni¬ 
que.  The  latter  method  has  been  used  to  check  the  load  distributions 
originally  obtained  by  the  former  so  that  an  independent  check  of  these 


results  has  been  carried  out. 

It  was  shown  that  the  lifting- 
surface  analog  to  the  solution  for 
load  distribution  over  an  unsteady, 
two-dimensional  wing  traveling  at  a 
constant  subsonic  speed  involved  the 
calculation  of  load  distribution 
over  a  swept— forward  wing  tip  with 
subsonic  edges.  Sketch  (ac)  indi¬ 
cates  the  geometry  associated  with 
the  boundary  conditions,  and  the 
solutions  are  calculated  for  the 
various  regions  shown. 


—I 


In  the  notation  of  the  unsteady 

problem  the  expression  for  the  velocity  potential  can  be  written 


cp  =  _i  rp  (vu)  ati  fa!  (A1) 

*  J  JT  J  (t-ti)2-(x-xx)2 

where  t  is  the  area  on  the  wing  plan  form  included  in  the  Mach  f ore- 
cone  from  the  point  (t,x).  Equation  (Al)  is  applicable' only  for  cases 
in  which  wu  is  known  at  all  points  within  the  forecone,  as  is  the  case 
when  the  edge  of  the  wing  within  the  forecone  is  everywhere  supersonic 
(i.e.,  region  1  in  sketch  (ac)).  However,  E ward  (reference  8)  has 
extended  the  solution  provided  by  equation  (Al)  to  include  cases  such  as 
shown  in  sketches  (ad)  and  (ae)  in  which  the  forecone  intersects  a  sub¬ 
sonic  edge  and  includes  a  region  of  unknown  upwash.  As  was  pointed  out 
in  reference  8,  equation  (Al)  applies  in  these  instances  if  the  area  of 
integration  t  is  limited  to  the  shaded  regions  shown  in  the  sketches. 
It  is  apparent,  therefore,  that  the  potential  (and  thus  the  loading) 
over  a  sinking  or  pitching  wing  can  readily  be  determined  for  regions  1, 
2,  and  k  in  sketch  (ac). 
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Points  in  regions  3  and  5  have  forecones  which  intersect  two  sub¬ 
sonic  edges,  and  the  method  just  discussed  can  no  longer  be  directly 
applied.  In  reference  8,  however,  a  method  was  given  of  evaluating  the 
upwash  in  the  region  between  the  Mach  cone  from  the  apex  and  the  leading 
edge  (region  6  in  sketch  (af)).  Thus,  the  plan  form  has  become,  effec¬ 
tively,  one  such  as  shown  in  sketch  (ag)  in' which  only  one  edge  is  sub¬ 
sonic.  This  reduces  the  problem  of  finding  the  potential  in  these 
regions  to  the  same  problem  as  was  involved  in  region  4.  The  analysis 
used  in  finding  the  loading  over  the  various  regions  will  now  be 
considered. 


(af)  (eg) 

First,  introduce  a  new  coordinate  system  in  which  the  lines  x  =  — t 
and  x  =  t  are  taken  as  the  r  and  s  axes,  respectively,  (see 
sketch  (ah)).  (This  amounts  to  a  rotation  of  the  original  axial  system 
through  an  angle  of  45°.)  The  transformations  relating  the  r,s  to  the 
x,y  system  are 
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r  =  (t-x)  t  =  -h  (r+s) 

a/2  JIT 


s  =  -pzz  (t+x)  x  =  -4r  (s-r) 

JT  J  2 


In  the  new  coordinate  system  equa¬ 
tion  (Al)  is  written 


(wu)dridsi 

'(r-rx)(s-s1) 


The  vertical  induced  velocity  • 
wu  over  the  wing  plan  form  is  given 
in  equations  (5)  and  (6)  for  the  sink¬ 
ing  wing  and  pitching  wing,  respec¬ 
tively.  The  method  developed  in 
reference  8  was  used  to  obtain  the 


(ah) 


value  of  wu  over  the  area  between  the  lines  s=0  and  r=s(l4l40)/(l--M0) 
(region  6  in  sketch  (af).  The  results  for  the  sinking  and  pitching  wing 


are,  respectively, 


In  terms  of  the  r,s  coordinate  system,  the  expression  for  the 
loading  can  be  written 


AP  =  4  1  /d_  _d_  \ 

q0  -  /rW  bsj  9 


(A4) 


If  the  subsonic  trailing  edge  is  not  included  in  the  forecone  from  the 
point  at  which  the  loading  is  being  determined,  the  expression  for  cp 
may  be  substituted  directly  into  equation  (A4).  However,  if  the  trailing 
edge  is  included  in  the  forecone,  and  if  the  Kutta  condition  is  to  apply 
along  such  an  edge,  it  can  be  shown  that  the  equation  for  the  loading 
coefficient  assumes  the  form 


H 1 d  Hid  Hi  d 
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AP  =  _  2  P  dri  /  _S_  _S 
Q-o  itVoMo  J7  J  r_r  \  Sri  Se 


(Vu)dS; 

y  S-Si 


Sinking  Wing 


The  preceding  method  can  he  applied  to  the  sinking  wing  to  obtain 
the  following  integral  relationship  for  the  loading  over  the  various 

regions.  The  subscripts  on  ^  ^  ^  indicate  the  region  for  which  the 

particular  equation  applies. 


/  Ap\  _  2 

V  40 )x  ~ 


(s?  A)./ 


S  r-ri  S-S 


A  +  A 

Sr  Ss 


r(l-Mo) 

(1+Mo) 


dsx  rr  dri 
s-Si  J-sx  J  r-rl 


(  Ae 

\  1o, 


_a  /  rr  dri  f  _a  +  A  ^  rs 

nM°  ^-'Jo  </ r-ri  \  ^ri  /  J n(l-Mo) 

(1+Mq) 


r  dri  / A.  +  A 

'o  V  r-ri  \Sri  Ss 


ri(l-Mp) 

(1+Mo) 


f  dsx 


(ri— s1)-M0(n+s1 


—  arc  tan 


-r,  ^  ^ri  )J- 


rd  vr-ri 


•ri  V  s-sx 
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=  j-  r  dri  (—  y~\  r  ds 

^g0//4  «Mo  ^  J r_ri  V^n  Ss  )  J-Vl  J  S-SX 


a 


!(**) 
a  \  W 


dri  f  d  d 


nMoyr d  f^\*n  ^/4(1- Mo)  <✓ 

(1+Mq) 


dsi 


s~Sl 


f  (J-+*-')  f 

Jr,  J  r-ri  \  drx  Ss  /  J 


ri(l-Mp) 

^T-2asi 

JC 


vhere 


2s 


(ri-siJ-MoCr^j.) 


arc  tan 


2s 


(n-si)-Mo(ri+Sx) 


S — Si 


rj  = 


s(l+M0)-cyr2 

1-Mq 


Most  of  these  integrals  can  he  readily  evaluated  to  give 


1  (  £P 

a-  \  qQ 


k_ 

Mo 


(A6a) 


!/■  AP 
a  v  q 


_8_  / 

jtMo 


[  J«£2_  /,,1.P-X0  +  arc  tan  /  Mbto+xp 

\  1+Mq  V  Mot0+x0  v  t0-x0 


/ Mpto+xp  \ 
v  Iq— xo  / 


(A6b) 


i  =  _§_  /Tv^  + JL  , 

a  \q°4  Jt(l+Mo)  v  Mot0+x0  hMq  L 


2X0-^0(1-440) 

arc  sm  — - -  + 

(1+Mo)  tG 


arc  sin  £<l-*o)-to(l+Mo)  ] 
to(l-Mo)  I 


i  ^ 

a  \40 


=  -r-  arc  sin  /  -°  1+Mq1o 


to(Mo-l) 


(A6c) 

(A6d) 


38 


MCA  TN  2403 


MCA  TN  2403 


39 


1 

9  V  90/2 


8 


3«Mo 


3/2 


(Wo)  +  3 / ( tQ-X0 ) ( Mot0+X0 )  * 

<1+Mo>  VMoV^ 


1  /  Ap\ 

9  \W3 


1  ( Ap\ 

9  \90/4 


3(M0t0+x0)  arc  tan  /^to+-9- 

v  ‘to-Xo 


(ATb) 


=  4  I 


1  Mo(l-Mo)  /(to-x0) 


y [ l-( Mpt 0+x0 ) ] [ ( t0+x0 )  -  1]  + 


(Mptp+x0) 


arc 


tan  A  (MototEal  _  arc  tan  / 

y  (t0+x0)-i  y 


t0-x0 


Mot0+x0 


(A7c) 


=  JL  arc  tan 

itMo  \  °  °  °  y  t0+x0-i 


y  ( ■fco+XQ-l ) [ l-( Mot  o+x0 ) ] 

[  yi=^w  + 


(A7d) 


/Ap\ 

_  _8_ 

f  140  r 

\90y5 

it  Mq 

lyi^o2  L 

( M0t0+x0 )  tanh  1  Jl-(Mot0+x0) 


tQ  x0)(to+x0) 
2 


y ( t Q— Xq )  ( Mpt  p+Xp)  +  ^.^ryg 


/(tp-Xo)" 

^0^0+X0 


( Mptp+Xp )  arc  tan  J  ^3 


+Xp 

Ko 


arc  tan 


/Mpto+Xp 

y  t0-x0 
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Mp  -l  / (tQ— x0)(l— Mp)  1  1 

V  (t0+x0)(l+Mo)  J  J  c 


(A7e) 


where 


rdU-Mp) 


-2Mq  */r+s  i  - 


— 4si(r+Si)  +  (r^+s1)[r(l-M0)  +  s^l+Mo)  +  2si(1-Mq)]  «  ~| 

arc  tan . . .  ■■  . . . .  +  2  ~ 

2  v/,2s1(r+s1)(r-rd)[rd(l-M0)-s1(l+M0)]  -I 


4Mq  J r— rd  arc  tan 


rd(  1-Mq  )— s  i  (  1+Mq  )  j— 


+  _2Mo2_  (yr(i_M0)_Sl(i+M0)  + 


£Ji^r  (2^  +  6>/^-|_^arctan 

3  1+Mq  |_  s-sx  J  y—  7  (l-M0)(r-rd)  J 


LIFT  AND  PITCHING^OMENT  C0EFFIC3ENTB 


The  lift  and  pitching— moment  coefficients  may  be  obtained  by  suit¬ 
able  integrations  of  equations  (A6)  and  (A7)  and  are  given  in  the  time 
intervals  indicated  in  sketch  (aa)  by  the  following  expressions: 

Sinking  wing 


0  <  t0  <  -i— 

1+Mq 


(A8a) 


’  =  -  sfc  [  i  -  t  (1"Mo)  +  ^r2- (Mo_2) 


(A9a) 


<  tn  <  -L- 

1+Mq  “  -  1-Mo 
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^3t0(l-Mo2) 

1+Mo 


4l 


where 

C]%.f  = 


-a 


itMo 


arc 


tan  + 

J  2t<,(l+Mo)-2 


It  (  tn  — 


1+Mq 

4— 2t0(l+Mo) 


)  +  (1+Jh)s  -s/[2t0(l+Mo)-2][2  -  to(l-Mo2)]  + 


1+Mo 


arc  sm 


t0(l+Mo)-l 

t0(l+Mo) 


^od«o)-l  +  [2-to<l*Mo)J  tan 

nl-Moto  /  Ap  \ 

J- 2__to  U°i  ° 

1+Mo  0 

—  j  is  given  by  equation  (A6e) 


(A8b) 


8 

ttMo 


tpg(  18Mq+9Mq2)  2to(Mo-l)  2 

IS  1+Mq  (l+Mo)2 


arc 


{[ 

tan  i  Z-toil-Mo2)  +  r-^t02(l-Mo)g  +  t0(Mog-4Mo+3)  +  Mpg+2Mo-3  -| 
011  V2[t0(  1+Mq  )-l]  l  16  +  2(l+Mo)  +  2(1+Mq)2 


arc  tan 


1  +  [2-t0(  1— Mpg ) ]  [  tp(  1+Mq PlJ 

#/t0(l+Mo)-l  1+Mo 


to(l-6Mo+9Mog) 
16 (1+Mq) 


1+9Mq 
4(1+Mo)2 J 


L-Mot0 


-  t, 


+  yt0(i+Mo)-i 

1+Mq 

(x0+Mot0)(^^^)  dx0 


5t0(l-Mo) 

16 


5+ 3Mp 
8(  1+Mq  ) 


(A9b) 


1+Mq 

Pitching  wing 
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c. 

< 

1+Mo 

c 


where 

ciuj* 


^  =-^{to3[|  (1-Mo)3  +“Mo]  +  ^  to^l-Mo)2  -  i  to(l-Mo)  +lj 

(Alla) 


t0  < 


1-Mo 


8 


h'  =  Vto+toMo-1 


t0(-^M02+2MQ-3)  +  ^Mq2+TMq+3 

12(l+Mo)  3(1+Mq)2J 


/TO?  ( |  -  + 


g(-Mog+2Mo+l)  ,  (1-Mo) 


-  t, 


+  1 


tan  / - -  falL  +  - 

•J  t0+t0Mo— 1  \  ^  2 


arc 


arc  tan 


/tp+l-' fcoMp 

to+t0M0“l  _ 


—  tanh  1 

2  */  1-Mo? 


1-MpV  to+toMp— l\ 

k l+Mp/  \  t  o+l— t  0Mq 


)]}*/': 


1— MqIq 


Gg 

2  c 
1+Mo 


dx0 
(AlOb) 


G2  is  defined  under  equation  (A7e). 
1  _ Mo 


8 

JtMo 


tanh 


_a  //1-Mq\  /t0(l+Mo)-l\  _ 
\1+Mq/  \t0(l-Mo)+l/ 


arc  tan 


3  J  1-Mo2 

(Mo+1)(-Mo2+4Mo+1)  -  ~  (1-Mo)2  +  Y  (1-Mo)  “  |j 

/  1  _  /Mot03  .  1 

</  t0(l+Mo)-l  v  6 


t0(l-Mo)+l 


+  -  larctoytKH  + 


yt0(l+Mo)-l  (-^[(l-Mof^Mo]  +  t0(9Moy^-ll) 

°  *  l  24 ( 1+Mq)  72(1+Mq) 

+  yt02-(l-Mot0)2 


39Mo2+54Mq+19 

36(1+Mo)2 

1— M^t 


t02(2+Mo2)  +  t0Mo  -  8 

- 35 - 


t/-i+  *  y 


}  + 

(Allb) 
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APPENDIX  B 

DETERMINATION  OF  SUPERSONIC,  INDICIA!,  SECTION  LIFT  AND 
PITCHING  MOMENT  CURVES 
THE  LOAD  DISTRIBUTION 


In  the  case  of  the  unsteady  supersonic  wing  the  expression  for  the 
velocity  potential  may  he  readily  obtained  by  placing  the  values  of  vu 
given  by  equations  (5)  and  (6)  in  the  equation 


<P  = 


-==^=±=== 

</(t-t;L)2-(x-x1)2 


(Bl) 


where  t  is  the  area  on  the  plan  form  included  in  the  Mach  fore cone. 
The  loading  may  then  be  calculated  from  the  relationship  given  in 
equation  (4). 


Sinking  Wing 


The  load  distribution  over  the  regions  A,  B,  and  C  shown  in 
sketch  (aa)  are  given  by  the  following  expressions: 


l/AP  \  =  _4_ 

aK<lo/A  Mo 


(B2a) 


9IH 
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*&)  _ _ ii_ 

yMo2-i  Ln 


arc  cos 


SggSa  +  arc  cos  ('_&')]  (B2b) 

x+Mot0  aMq  \  t0  /  J 


i  = 


-  i‘ 

a  Vi 


q%  V  Mo2- 1 


(B2c) 


Pitching  Wing 


For  the  case  of  the  pitching  wing  the  values  of  ~  in  regions 

regions  A,  B,  and  C  are,  respectively. 


q  (  S)A  Mo  (xo+M°to) 


(B3a) 


1  f  Ap  \  __  4  Xp+Mpto 

q  \  q°  /jg  *  L  ^Mq2-! 


Mpx0+tp  ,  Xp+Mptp 
arc  cos  -~-u — ~  arc  cos 

*o+Mot0  Mq 


(-£) 


y  to2-3 


1/  ap\ 

q  ^  q  JL 


J  Mq2-1 


(xo+Moto) 


(B3b) 


(B3c) 


LIFT  AND  PITCHING-MOMENT  COEFFICIENTS 


The  lift  and  pitching-moment  coefficients  may  he  obtained  by  suitable 
integrations  of  equations  (B2)  and  (B3)  and  are  given  in  the  time 
intervals  indicated  in  sketch  ( aa )  by  the  following  expressions: 

Sinking  wing 


0 


<  t. 


<-U_ 

-  1+Mo 


cl 


a 


_4 

Mq 


(B4a) 


c 


(B5a) 


MCA  TN  2403 


45 


46 


MCA  TN  2403 
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TABLE  I 

TABULAR  VALUES  OF  SUBSONIC  INDICIAL 
LIFT  AND  PITCHING  MOMENT  CURVES 
M0  =  0.8 


to 

|3 

5c*a 

2(3  , 

3^8 

2p  » 

TCma 

2(3  . 

—  cm  * 
it  111  a 

0 

0.478 

0.318 

-0.955 

-0.637 

.1 

.468 

•314 

-.931 

-.618 

.2 

.458 

.312 

-.898 

-.601 

•  3 

.449 

•  313 

-.857 

-.587 

.4 

•  439 

•  317 

-.805 

-.578 

•  5 

.430 

.325 

-.745 

-•575 

.6 

.423 

•  336 

-.689 
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.8 

.428 

•  359 

-.639 

-.600 

1.0 

.442 

.383 

-.621 

-.620 

1.5 

•  479 

.438 

“  .608 

-.664 

2.0 

•  515 

.484 

-.613 

-.694 

2.5 

.542 

.525 

-.619 

-.716 

3.0 

•  574 

•  558 

-.625 

-.727 

3-5 

•  599 

.586 

-.635 

-•735 

4.0 

.619 

.610 

-.645 

-.740 

4.5 

•  637 

.630 

-.658 

-.747 

5-0 

.652 

.645 

-.671 

-.755 

6.0 

.678 

.672 

-  .696 

-•773 

7-0 

•  701 

.696 

-.718 

-.790 

8.0 

.722 

•  717 

-.737 

-.804 

9 

.740 

.736 

-.754 

-.817 

10 

.757 

•  752 

-.769 

-.828 

15 

.818 

.815 

-.827 

-.871 

20 

.858 

.855 

-.865 

-.899 

25 

.885 

.883 

-.890 

-.918 

30 

.904 

•  903 

-.908 

-•932 

40 

.929 

.928 

-•932 

-.949 

50 

.945 

.944 

-.947 

-.961 

6o 

•  955 

•  955 

-.957 

-.968 

80 

.968 

.968 

-.969 

-.977 

100 

•  976 

•975 

-.976 

-.982 

CO 

1.000 

1.000 

-1.000 

-1.000 

VQt /c  t  Chord  lengths  traveled 
(o)  Lift  on  a  sinking  wing. 


Figure  /.-  Variation  of  two-dimensional  indicia /  lift  and 
pitching-moment  coefficients  with  chord  lengths 
traveled  for  several  Mach  numbers. 
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0  2  4  6  8  10 

V0t'/c  ,  Chord  lengths  traveled 


(b)  Pitching  moment  (about  leading 
edge)  on  a  sinking  wing. 


Figure  l.  -  Continued. 
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(c)  Lift  on  a  wing  pitching  about  its  leading  edge. 


Figure  /.  -  Continued. 
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Vj'/c ,  Chord  lengths  traveled 

(d) Pitching  moment  (about  leading  edge)  on  a 
wing  pitching  about  its  leading  edge. 


Figure  I.-  Concluded. 
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(a)  Sinking  wing. 

Figure  2-  Variation  of  two-dimensional  load  distribution 

with  percent  chord  for  a  Mach  number  equal  to 
0.0. 
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Percent  chord 

(b)  Wing  pitching  about  its  leading  edge. 

n 


Figure  2.-  Concluded. 


(a)  Sinking  wing. 

Figure  3. -Variation  of  two  -  dimensional  indicia /  load 
distribution  with  percent  chord  for  a  Mach  num¬ 
ber  equal  to  t.2 . 


<5* 


56 


MCA  TN  2403 


Percent  chord 


(b)  Wing  pitching  about  its  leading  edge. 
Figure  3.  -  Concluded. 
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